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New types of stationary solutions of a one-dimensional driven sixth- 
order Cahn-Hilliard type equation that arises as a model for epitaxially 
growing nano-structures such as quantum dots, are derived by an exten- 
sion of the method of matched asymptotic expansions that retains ex- 
ponentially small terms. This method yields analytical expressions for 
far-field behavior as well as the widths of the humps of these spatially 
non-monotone solutions in the limit of small driving force strength which 
is the deposition rate in case of epitaxial growth. These solutions ex- 
tend the family of the monotone kink and antikink solutions. The hump 
spacing is related to solutions of the Lambert W function. 

Using phase space analysis for the corresponding fifth-order dynamical 
system, we use a numerical technique that enables the efficient and accu- 
rate tracking of the solution branches, where the asymptotic solutions are 
used as initial input. 

Additionally, our approach is first demonstrated for the related but 
simpler driven fourth-order Cahn-Hilliard equation, also known as the 
convective Cahn-Hilliard equation. 

Keywords: convective Cahn-Hilliard, quantum dots, exponential asymptotics, 
matching, dynamical systems 

1 Introduction 

A paradigm for phase separating systems such as binary alloys is the Cahn- 
Hilliard equation for the phase fraction u 
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where Q{u) is the derivative of the double- well potential J 7 , typically 



Q(u) = T'(u) =u — u. (1.2) 

The long-time dynamics are characterized by the logarithmically slow coarsen- 
ing process of phases, corresponding to local minima of the potential, separated 
by interfaces of width e. This process is well described by the motion of equidis- 
tantly spaced smoothed shock solutions or kinks ( "positive kinks" ) and antikinks 
("negative kinks") which connect the local minimum of T(u) at it = —1 to that 
at u = 1 and vice versa. 

In recent years, an extension of this model has been studied for the case 
when the phase separating system is driven by an external field 1(1, 27 1 . In one 
space dimension it can be written as 

Ut - vuu x + (Q(u) + s 2 u xx ) xx = 0, (1.3) 

where v denotes the strength of the external field. This equation, the convec- 
tive Cahn-Hilliard (CCH) equation, also arises as a model for the evolution of 



the morphology of steps on crystal surfaces 2l|, and the growth of thermody- 
namically unstable crystal surfaces into a melt with kinetic undercooling and 
strongly anisotropic surface tension @,[ll, 17 [. 



The dynamics of this model as v — > are characterized by coarsening, as 
is typical for the Cahn-Hilliard equation {y = 0) 0, 0]. If v — > oo using 
the transformation u h- > ujv in (|1.3|) one obtains the Kuramoto-Sivashinski 
equation, which is a well-known model for spatio-temporal chaotic dynamics 
(see e.g. [lj| and references therein). Recently, Eden and Kalantarov Q also 
established the existence of a finite-dimensional inertial manifold for the CCH 
equation, viewed as an infinite-dimensional dynamical system. 

A related higher-order evolution equation arises in the context of epitaxially 
growing thin films (for a review on self-ordered nano-structures on crystal sur- 
faces see Shchukin and Bimberg j24[). Here, the formation of quantum dots and 
their faceting has been described by the sixth-order equation 

Ut - vuu x - (Q(u) + £ 2 u xx ) xxxx = 0, (1.4) 

where u denotes the surface slope, v is proportional to deposition rate [22| 
and Q(u) is given with (|1.2j) . it is assumed to have this form from now on 
throughout the paper. The high order derivatives are a result of the additional 
regularization energy which is required to form an edge between two plane 
surfaces with different orientations. This implies that the crystal surface tension 
also depends on curvature, which becomes very high at edges as the parameter e 
goes to zero. In analogy to the Cahn-Hilliard equation, here the phases are the 
orientations of the facets. This higher- order convective Cahn-Hilliard (HCCH) 
equation shares many properties with the CCH equation. In both cases the 
dynamics are described by conserved order parameters if v = 0. They also 
share characteristic coarsening dynamics as v — * and chaotic dynamics as 
v gets large. To understand the complicated structure of the solutions it is 
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instructive to study first the stationary solutions and their stability as it has 
been done for the CCH equation 2, 231 ■ For small v the stationary solutions 



for both equations have been characterized by the monotone kink and antikink 
solutions \wL |22| . Recently new spatially non- monotone solutions were found 
for the lower order equation [28[. In this study we establish that the HCCH 
equation also possesses such non-monotone solutions. We show this by using 
phase-space methods for the corresponding fifth-order boundary value problem. 
We use the expression "simple" or "monotone" for a solution that connects 
the maximal value of u(x) to the minimal value without any humps on the 
way down, although these extrema exist and lead to non-monotonicity even for 
simple (anti-)kink solutions of the HCCH equation. 

Since the treatment of this high-order problem is not straightforward, one 
part of this study is concerned with the development of an approach that accu- 
rately locates the heteroclinic connections in the five-dimensional phase space. 
We find that these stationary solutions develop oscillations whose width and 
amplitude increase as v — > 0. 

In the second part of this study we derive an analytic expression for the width 
and amplitude within the asymptotic regime of small external field strength. For 
the CCH equation we find that the width has a logarithmic dependency on the 
strength of the external field, while for the HCCH equation our analysis yields 
a dependency in terms of the Lambert W function. In order to arrive at these 
expressions we solve the fifth-order equation by a combination of the method of 
matched asymptotic expansions and exponential asymptotics. We first demon- 
strate our approach for the third-order boundary value problem arising from 
the CCH equation. Our approach generalizes the work by Lange [IH to higher- 
order singularly perturbed nonlinear boundary-value problems, where standard 
application of matched asymptotics is not able to locate the position of interior 
layers that delimit the oscillations of the non- monotone solutions. 

Reyna and Ward (lij previously developed an approach to resolve the in- 
ternal layer structure of the solutions to the boundary-value problem for the 
related Cahn-Hilliard and viscous Cahn-Hilliard equations. The approach is 
based on a method due to Ward [25| who uses a "near" solvability condition for 
the corresponding linearized problem in his asymptotic analysis, and who was 
inspired by an earlier variational method [l3| and work by O'Malley [IH and 
Rosenblat et al. , who investigated the problem of spurious solutions to singu- 
lar perturbation problems of second-order nonlinear boundary-value problems 

Moreover, for the related Kuramoto-Sivashinsky equation, a multiple-scales 
analysis of the corresponding third-order nonlinear boundary- value problem by 
Adams et al. [l| shows that the derivation of monotone and oscillating traveling- 
wave solutions involve exponentially small terms; their method is based on an 
analysis of the Stokes phenomenon of the corresponding problem in the complex 
plane (see Howls et al. [l2[ for an introduction). 

In what follows we begin with the phase space analysis for the CCH equation 
in section 2, followed by the asymptotic treatment for v« 1. The asymptotic 
ideas used for the CCH equation are then applied to the HCCH equation in 
section 3. The solutions obtained there are useful to serve as initial input for the 
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numerical investigations of the branches of non-monotone solutions in section 4. 
In this part we develop our numerical approach and then use it to identify new 
stationary solutions of the HCCH equation, these agree with the asymptotic 
theory. Finally we shortly sum the results up together with concluding remarks 
in the last section 5. 



2 Stationary solutions of the convective Cahn- 
Hilliard equation 

The high-order term in the CCH equation represents the regularization of the 
internal layers of the solutions. For most of our investigations we consider 
the problem in the scaling of the internal layers, or inner scaling, where the 
i-coordinate is stretched about the location x — x of a layer according to 

(2.1) 



e 

In this scaling the CCH equation becomes (after dropping the "*") 

e 2 u t - -{u 2 ) x + (Q{u) + u xx ) xx = 0, (2.2) 

where S = ev. The stationary problem obtained by setting u t to zero can 
be integrated once, requiring that the solutions approach the constants ±vG4 
as x — > ^foo, where A is a constant of integration. That is, we consider the 
boundary value problem 

- (u 2 - A) = (Q(u) + u xx ) x (2.3) 

together with the far-field conditions 

lim u = (2.4) 

x — >±oo 

and vanishing derivatives in the same spatial limit. We refer to solutions of 
this system as antikinks. Monotone antikinks are known analytically (l6j . while 
recently, non-monotone connections were computed numerically by Zaks et al. 
[jjj ]. We now briefly discuss the numerical approach to obtain these solutions. 
Here we concentrate on the regime where < 5 <C 1 in order to compare with 
the asymptotic solutions derived later on. For a bifurcation analysis for larger 
S we refer to [281 ] . 



2.1 Numerical solutions 

For the numerical solutions we will work with a rescaled version, where we set 
u = \f~A~c so that the equilibrium points do not depend on A, and for Q(u) given 
by (fT2]h (j2~3l) and become 

2 

(1 - c 2 ) = -—={c xx + c- Ac 3 ) x , lim c=Tl- (2.5) 
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For this problem we find it most convenient to present a shooting method which 
enables us to track solution branches in the (A, 8) parameter plane. We trans- 
form l|2.5[) to a first order system U' = F(U), where F : R 3 — > II 3 is the 
function 



F 1 (U)=U 2 , F 2 (U) = U 3 , F 3 (U) = (3A(U 1 ) 2 -l)U 2 + ^((U 1 f-l), 

(2.6) 

We work in a three dimensional phase space and denote cither vectors or whole 
trajectories therein with capital U's. We use the same notation for two different 
objects, because it will be clear from the context what is meant. Subscripts 
indicate the components. 

The characteristic polynomials at the equilibrium points U = (±f,0,0) T 

are 

A 3 + A(l - 3A) t <SV^ ■ (2.7) 



The signs of the real parts of the roots determine the dimension of the stable and 
unstable manifolds W u (U+) 7 W S (U~), W S (U+), W U (U~) of the equilibrium 
points. The latter two are two-dimensional and so the existence of a kink is 
generic, while this is not the case for the antikinks. The dimensions of W U (U + ) 
and W S (U~) are one, so that the heteroclinic connections from the positive 
to the negative equilibrium arise from a codimension two intersection. This 
means that with the two parameters A and 8 we can expect only separated 
solutions when the manifolds intersect, but due to the reversibility properties 
which are discussed below the codimension reduces to one and we can expect 
separated solutions for the free parameter A and a fixed 8, hence one or several 
whole branches in the (A, 8) parameter-plane. An example of a non-monotone 
connection is sketched in figure [T] where the trajectories wind themselves in the 
phase space with a solution that exhibits f5 humps. 



Reversibility and computations It is instructive to note that the solution 
of (|2.5|) is translation invariant, c(x) — > c(x + L), and forms a reversible dynami- 
cal system, hence the solutions are invariant with respect to the transformation 
x — > —x, c — > — c, as has also been noted by Zaks et al. pjj . 

Let us consider generally a fc-dimcnsional phase space, since the following 
discussion will be also useful in section [4] where we analyze the HCCH equation 
with its higher order system. The linear transformation 

R : M k -> R fc , R{U.j) = (-l) j Uj,j = 1, . . . , k (2.8) 

fulfills R 2 = Id and RF(U) = -F(RU) for k = 3 and (J2U) and represents 
the reversibility in the phase space. It is an involution (or a reflection) and its 
set of fixed points is the symmetric section of the reversibility, these are zero 
at odd components, Ui = for odd i. A solution that crosses such a point 
necessarily symmetric under R, and for each point U on the connection there 
exists a corresponding transformed point RU somewhere on the branch. In 
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Figure 1: CCH: Antikink solutions connecting the hyperbolic equilibrium points 
U + and U~ are sought in a 3-D phase space. The unstable manifold emerging 
from U + , W U (U + ) is one-dimensional, as is the stable manifold W S (U~). The 
approximating linearized spaces E U (U + ) and E S (U~) are drawn as dash-dotted 
lines and are used in the computations. 

fact there is an equivalence here since odd solutions necessarily cross a point in 
the symmetric section. It holds that c and its even derivatives have to vanish 
in the point of symmetry L because of the fulfilled equations £ im c(x + L) = 

— j Tm. c{—x + L),m = 0,1,..., L^/2J, and continuity of the solution and its 
derivatives. 

From the above we conclude that with a shooting method we can stop inte- 
grating when we find a point with zero odd components, since the second half 
of the solution is then given by the set of transformed points under R. Hence 
we define the following distance function for a trajectory U over the interval of 
integration which helps to find these points 




Figure 2: Distance function cLa defined by (|2.9|) depending on A with fixed 
5 = 0.05, showing the first 14 zeros corresponding to heto to hetn. 

The minimization of dj\(U) over the free parameter. mind^fE/), must result 

A 

in the value zero for an anti-symmetric heteroclinic solution. We can use this 



G 




-11 -8 -5 -40 40 

log (8) x 



Figure 3: Parameter plane, log (5) for the x- and log(l — A) for the y-axis, for the 
CCH equation for the first 9 antikink solutions hetk, k = 0, 1 . . . , 8. The graphs 
on the right show the shapes of representative het^ solutions, hence those on 
the fifth line from below, for the approximate (A, 5) tuples (0.8259,0.0289), 
(0.9893,0.0017) and (0.9998,2.6457- 10~ 5 ). 

condition for shooting and BVP formulations, for both the CCH and later the 
HCCH equation in section 2] 

For a fixed value of 6 and a range of different A we follow the relevant branch of 
W U (U + ) by shooting from an initial point U + ± ev near U + , where v is a unit 
eigenvector corresponding to the positive eigenvalue of dF/dU\ u=u + and e<l. 
We stop the integration if a certain threshold value for \Ui \ is crossed. Figure [2] 
shows cIa(U) as a function of A for S = 0.05. 

At this point we have heteroclinic connections for one fixed value of 6 which 
we denote by hetk,k = 0, 1, . . . (using the notation in Zaks et al. [H|). heto 
is the analytical, monotone tanh solution while hetk has k humps on the way 
down from \f~A to — \f~A. We will use the same terminology for the solution 
structure of the stationary HCCH problem in section [4] Here, a hetk solution 
corresponds to the fcth zero from the right in figure O We then follow the 
roots of the distance function by linearly extrapolating to a new guess for A 
and use a bisection algorithm to converge fast to the next root. Figure shows 
a portion of the (A, 5) parameter-plane, where we concentrate on very small 
values of <5, or differently interpreted, on the bifurcation of the various spiraling 
CCH orbits from the heteroclinic connections of the CH equation in its one 
dimension smaller phase space. 
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2.2 Asymptotic internal layer analysis 

For the asymptotic analysis we use a slightly different scaling than for the nu- 
merical treatment. Here, we let 



(2.10) 



denote the inner variable about a layer located at x — x. For the stationary 
problem we then obtain 



(u" + 2Q(u))' = SV2(u 2 - A) 



(2.11) 



instead of (|2.3p . where ' = d/dx* . For later comparisons of the numerical and 
asymptotic results we have to keep in mind that the spatial scales differ by a 
factor of V2- 

We point out that the problem considered here shares the internal layer 
structure of the singular perturbation problems discussed by Lange [l5j |. and 
we will make use and extend this ansatz for our situation. This will also prove 
useful to understand the approach taken for the HCCH problem in section 13. li 
since there we have to carefully combine the exponential matching with the 
conventional matching procedure when matching the two regions. For both 
problems the asymptotic analysis can be conveniently carried out in terms of 
the small parameter 5. 

In the following analysis we consider the simplest case of a non-monotone so- 
lution with only one hump, as illustrated in figure[4l we note that non-monotone 
solutions with more oscillations can be treated similarly. 



2.2.1 The 1-hump solution 



outer region 
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Figure 4: Sketch of a 1-hump, or het\ solution showing the general setup for 
the matching procedure for the CCH and HCCH equations. 



We observe that the 1-hump solution has three internal layers, one at n m < 0, 
one at n p and one at the symmetry point in between. Since the solution is point 
symmetric we can choose this point to be x = and it will be enough to only 
discuss the two layers at n m and zero and then match them to the outer solution. 



Internal layer near n rn For the first internal layer at n m we let 
where R m < and set 

oo 

Km = K m + V2^ S k K mk , (2-13) 
fe=l 

so that to leading order the location where the solution crosses zero is R m and the 
additional terms account for the corrections due to the higher order problems. 
With u, n (x m ) — u(e(R m + s/2x m )) the governing equation becomes 

uZ + 2Q'(u m ) =8V2(u m - A), where / = (2.14) 

For the boundary condition where u m crosses zero we have 

'^^)=0 (2.15) 

and the condition towards — oo is 

lim u m (x m ) = VA. (2.16) 

x 1Tl — oo 

We now assume u m (x m ) can be written as the following asymptotic expan- 
sion, valid near n m 

oo 

U a (Xa) = U a o(x a ) + ^ 8 k U ak (x a ) , ( 2 -17) 
fe=l 

with a = m here. Additionally, we assume A has the asymptotic expansion 

A= l + SAi +0(5 2 ). (2.18) 
Observe that from (f2~16|) and |2~T8)) 

lim u m (x m ) = lim u m0 (x m ) + } 5 k u mk (x m ) = 1 + - S k A k . 

X m — > — OO Xjn — > — OO * * J, L 

k=l k=l 

(2.19) 

To leading order in 5 we get the problem for the Cahn-Hilliard equation 

uZ + 2Q'(u m0 ) = (2.20a) 
!i m o(0) = and lim u m0 (x m ) = 1 (2.20b) 

Xm — > — OO 
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with the unique solution u m o(x m ) = — tanh(x m ). Next, the problem of order S 



is 



(C(u ml ,x m )^ =V2 (tanh 2 (a; m )- 1) (2.21a) 
Umi(0) = Kmi and lim u m x{x m ) = — ^ , (2.21b) 

x m — > — oo Z 

where K m i and are constants to be exponentially matched and the operator 
£ is defined by 

C(v, z) = v" + 2 (1 - 3 tanh 2 (z)) v , (2.22) 

and z = x m , v = u m i and / = d/dx m . Note that the first boundary condition is 
obtained by expanding (|2.15p 

u m (^f^S k K mk ^j = u m (Sn ml +S 2 K m2 + 0{5 3 )^ (2.23) 

= Wmo(O) + S( K m iu' m0 (0) + Uml(O)) + 0{5 2 ) 

so that collecting the terms of order <5 gives 

"ml(O) = -K m x U m0 (0) = K m l ■ 

Next, we integrate (|2.21[) once to obtain 

C(u m x,x m ) = f m (x m ), (2.24) 

where f m (x m ) — — v / 2tanh(a; m ) +c m . Taking the limit of this equation to — oo 
yields c m = — v2 — 2Ai so that 

f m (x m ) = -\/2(tanh(a; ro ) + 1) - 2A 1 . (2.25) 

The homogeneous solutions of (|2.24p are 

<pm(x m ) = -u' m0 (x m ) = 1 - tanh 2 (.T ln ) , (2.26) 

ipm(x m ) = ( [ ,f Z f \ ] 4>m(x m ) ■ (2.27) 

Also note that linx^^-oo (j) m {x m ) = and ip m (0) = 0. At this stage it is conve- 
nient to choose the inhomogencous solution that remains bounded MS -I in ^ — OO 
and vanishes at x m — which is satisfied by 

(f> a f a dz - <p a (x a ) ip a f a dz (2.28) 

- oo JO 

with ol = 77i. Hence, the unique solution for (|2.2ip is the linear combination 

U m l(x m ) = -K m l4> m {x m ) + fm(Xm) ■ (2.29) 



10 



Internal layer near x = For the internal layer near the origin we proceed 
as above. Here, we stretch the independent variable as 

and construct an asymptotic expansion (|2.17p near x = with a = for the 
solution of the problem 

Uq + 2 Q'(u ) =SV2(ul- A), where / = (2.31) 

axo 

We note that the point x = is assumed to be the symmetry point of the 
complete solution, hence here we require 

it (0) = and <(0) = 0. (2.32) 

In anticipation of the exponential matching we also require that liniz^-oo uqo(xo) 

— 1, so that the solution to the leading order problem is uqo(xq) — tanh(xo) For 
the solution to 0(8) we find 

uqi{x q ) = b a ^ (x ) + ip (x ) (2.33) 

where bo is a further constant to be exponentially matched. Here, the homoge- 
neous solutions are 

<M^o) = -w'oo^o) and Vo^o) = \j ^T(7j J "M^o) (2.34) 

and the inhomogeneous solution is defined by (|2.28p . where a = and fo(%o) = 

— \/2tanh(xo). They are chosen such that (po(0) = and <£o(0) = 0, in fact wc 
have lim X0 _, ±oo ip {x ) = ±\/2/4. 

2.2.2 Exponential matching 

Exponential matching requires that all exponentially small and exponentially 
growing terms have to be accounted for and matched. This means first that 
we have to express the variable xq in terms of x m (or vice versa). From the 
definitions of these variables it follows that 

%o = x m + . (2.35) 

In particular, exponential terms in the solution uq(xq) transform as e 2x ° = 
e \/2K m £ 2x m an( j so f or ^jj f or n ig ncr order exponential terms e 2nx " or terms with 
different signs in the exponent. 

Now note that as xo —> — oo the leading and 0(S) solutions can be written 

as 

uoo(xo) = -1 + 2e 2xo - O{e 4xo ) (2.36) 
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and with \i = (|6q + V%) 



«oi(xo) = - ^ - ^e" 2 ^ + (H 6o + _L + e 2x + o(e 4, 0) _ (2 37) 
Written in x m variables the solution 

u (x m ) = -i + 2e 2x ™e^ R ™ + 0(e 2 ^ 2ft ") (2.38) 

4 M lg e + 46 0+ ^2 M m V2e Jj 

+0(e 2 ^"-) J + 0(,5 2 ) 

has to be exponentially matched to 

u m {x m ) = -l + 2e- 2x ™ + 0(e~ 4Xm ) (2.39) 
+5 (-(A, + ^) - i(A x + -±=)e 2 ^ + + ^V2 + 4« ml )e" 2 ^ 

-(3^i + w" 2 *'" + 0(e- 4 ^) ) + 0(<5 2 ) 

While we have already anticipated matching of the constants 
during the derivation of the leading order solutions, the constant terms of the 
0(S) solutions are first to be matched. Matching to the exponential terms in 
(|2.39[) entails a rearranging of terms of different orders of magnitude in the 
expansion (|2.38[) . In particular, the first exponential term to leading order 
in ()2.39|) matches the second term of 0(S) in (|2.38|) . the second and largest 
exponential term of 0(S) in (|2.39|) matches the second term of the leading order 
in ()2.38|) . and so forth. Summarizing, we obtain 

!(!„„ + v^) = .4, + -^ = 2, + 2. (2.40, 

where we denote p = 5 e~^ Km . Solving yields 

p = AV2, Ai=-4= and b = — . (2.41) 
v 2 v 2 

We observe that we have determined the 0(5) correction A\. Additionally, we 
now know that <5e~ v/2Km = 4^2, hence 

and if we recall (|2.13p and K m < then the width of the hump is —K m , where 

^ = ii|)_MV2) +0(i) (243) 
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Figure 5: Distances between the first two roots of the het\ solutions versus 
log (<5) together with the width predicted by the asymptotic formula (|2.43p . 

Further constants, such as K m i are found by including higher exponential terms 
and expansions of the higher order problems. Finally we note that making use of 
the symmetry of the solution about the point x = 0, the exponential matching 
of the solution near zero to the one near k p proceeds analogously. 

2.2.3 Comparison of numerical and asymptotic solution 

For the comparison with the asymptotic solution we are interested mainly in the 
heti solution which we derived in section [2.2.11 By numerical continuation of the 
shooting method, one obtains TV tuples (A^,S^),j = 1, . . . , TV in the param- 
eter plane that give a Tieii-branch when being connected. We use two vectors 

of parameters which we abbreviate A = {A^)j = \ t ... t N and 8 = (5^')j—i n 

to confirm the formulas we obtained in the previous section. Further we make 

use of a distance vector K = (K^)j-i jy K contains the distances between 

the zero crossings of the solutions, or in context of the asymptotics section (see 
figured} K*j) r* \K m {8^')\. To obtain the relation between A and 6 and the 
evolution of the distances we solve the least squares problems 

min ||(1 - fixS) - A\\l and min ||r?i log(<5r? 2 ) - , 

Mi »7i,»72 

hence we assume a linear law for the ^4-values in 6 and a general logarithmic 
law for the distances. We obtain 

A « 1 - 2.125 w 1 - -^=6 and K w -0.71 log(0.18<5) (2.44) 
v2 

which confirms the results from the analysis (|2.41j) and (|2.43|) . We see the good 
match in the distance plot in figure [5] 

These results motivated us to obtain a general rule for the relation between 
the two parameters of the CCH equation for different stationary solutions. The 
numerically computed branches in figure [3] show that the slopes of the hetk 
branches are one when plotting log(5) against log(l — A), so that the relation 
log(<5) + const = log(l — A) shows the linear dependence A(k) = 1 + Ai(k)S, 
where A(k) is the A value for the hetk solution and Ai{k) its linear coefficient. 
We see that the dependence of A\ with respect to the order of the heteroclinic 
connection k behaves linearly and we obtain a general expression for the squared 
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far field value of non-monotone hetk solutions, it is given with 

Al (k) = -^1. (2.45) 



3 Matched and exponential asymptotics for the 
stationary HCCH equation 

As for the CCH equation we will perform our analysis of the internal layers 
in the inner scaling (|2.10[) . From the stationary form of (|1.4p we obtain the 
equation 

(tt" + 2 Q(u))"' = -5 2 3 / 2 (u 2 - A) , (3.1) 

after integrating once and requiring that for an antikink \im x —,± 00 u = ^fy/A 
and setting 6 = e'v here. We consider the heb\ (one hump) solution, and again 
make use of the point symmetry of the problem. Now however, unlike for the 
CCH equation, the solutions in the outer region are not just constants. Here, we 
have to introduce an outer layer to the left of the inner layer about K m , see also 
figure 2] for the case of a 1-hump solution. In the following subsections we first 
briefly derive the solution to this outer problem and match it to the solution to 
the inner problem The remaining degrees of freedom arc then used to 

exponentially match it to a second inner layer near x = 0. 



It has been demonstrated in [22| for monotone antikink solutions of the 
HCCH equation, that it is necessary to match terms up to order 5 in order to 
obtain the correction A\, given the asymptotic expansion of A 

oo 

A=l + J2 5h/3 A k- (3-2) 

fc=i 

Here, for the non-monotone antikinks we have to match inner and outer solu- 
tions and then also exponentially match the inner layers. This has to be carried 
through iteratively up to three orders of magnitude in order to obtain not only 
the correction A\ but also the expression for the width of the humps. 

3.1 The 1-hump solution for the HCCH equation 

We start by shifting to the inner coordinates that describe the region near K m , 
which is to be matched to the outer region. Again defining x m by (|2.12[) . the 
governing equation in this inner region is 

(ul + 2Q(u m )) = - 2 3/2 5 (u 2 n - A) where /=— . (3.3) 
V / dx m 

For the boundary conditions we again place K m near the point where u m crosses 
zero, i.e. 

(3.4) 
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The condition towards — oo is not as trivial as for the CCH equation but needs 
to be matched to the outer solution in the region to the left of K m (or to the 
right of k p , taking account of symmetry) . 

For the outer region (see figured]), where x m becomes very large, we use the 
ansatz 

S = 5 1 ' 3 x m and Y(£;6)=u m (x m ;6) (3.5) 
and obtain the outer problem 

(> 3 Y U + 2 Q (Y) ) = - 2 3/2 {Y 2 - A) (3.6) 

with the far field condition 

lim Y{£) = VA. (3.7) 

f; — > — oo 

The region near x = 0, for which we use the variable Xq from (|2.30[) . is 
described by the problem 

(u^ + 2Q(u Q )) = - 2 3 / 2 S (uq - A) where /=— (3.8) 
V / dxo 

The point x — is the point of symmetry of the solution. Here we require 

u (0) = 0, u'o(0) = and u' '"(0) = , (3.9) 

plus additional conditions from the exponential matching to the internal layer 
near K m as xq — > — oo, as we have shown for the CCH equation. 

Here we assume the solutions to these three problems for Y, u m and uo can 
be represented by asymptotic expansions 

oo 

u a (x a ; e) = u a o(x a ) + <5 fe / 3 u a k(x a ), where a = 0,?7i (3.10) 

k=l 

valid near K m and x = 0, respectively, and 

oo 

Y(£; 8) = + E 6k/3 ( 3 - n ) 
fc=i 

valid in the outer region, where we let 

oo 

K m = K m + V2^2S k/3 K mk . (3.12) 

fe=l 

To obtain solutions to the outer problem is straightforward [22| , but in order 
to be more comprehensible we include the results in appendix [X] The solutions 
to the other regions are discussed now. 
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3.1.1 Leading order 

To leading order in S we get the problem 

(u^ + 2Q(« ro0 )) / " = (3.13a) 
u mQ (0) = (3.13b) 

Matching to the leading order outer solution (|A.2j) Yq = 1 we find 

UmO^ro) = - tanh(x m ) . (3-14) 

Its representation towards the internal layer about x = is given by 

u m o = -1 + 2e- 2x ™ - 2e- 4x ™ + 0{e~^) (3.15) 

The leading order problem for this region is 

«, + 2Q(«oo))'" = (3.16a) 

m 00 (0)=0, u' O q(0)=0 and it^'(O) = (3.16b) 

and its solution is 

Moo(^o) = tanh(xo) . (3-17) 
As Xq — > — oo its behavior is given by 

it 00 = -1 + 2e 2xo - 2e 4x ° + 0(e 6x °) . (3.18) 

3.1.2 0(S 1/3 ) 

Internal layer near x = K m The expansion of (|3.3[) and (|3.4[) to order J 1 / 3 
yields 

£(u m i,a; m ) = f m i(x m ) (3.19a) 
Umi(O) = -u^o(0) « m i = K ml (3.19b) 

where C is defined by (|2.22[) as for the CCH equation and 

/rnli^rn) ■ — Ci m X m -f- C2m^ra ~t~ ^3m - (3.20) 

The homogenous solutions are therefore (|2.26|) and (|2.27| . The constants Ci m , C2 m , 03^ 
are obtained by three successive integrations of the ODE for u m \ obtained at 
this order. We choose the inhomogeneous solution so that it grows only alge- 
braically and vanishes at x m = 0. Particular solutions to (|3.19b)) 
are of the form 



(P(Xj (X a 

) = 4la{x a ) / 4>a faj dz-(j} a (x a ) I 1p a f aj dz + J a j1p a (x a ) , (3.21) 

Jo Jo 
so that now we obtain tp m i for a = m, j = 1 in (|3 . 2 1 [) and 

7mi = -jrCim + ln(2)c 2m - c 3m . (3.22) 
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Hence the solution is 

u m i(x m ) = -K m i<f) m (x m ) + tp m i(x m ) . (3.23) 

We evaluate ip a , 4> a etc. and subsequent functions with the assistance of Maple. 
As the limiting behavior of u m i is 

u m i(x m ) = ~(ci m + 2c 3m ) - jc 2m a; m - ^ci m a;^ (3.24) 

tt(-7ci to - 8c 3m + 256k to i + 30c 2m + 4c 2m 7r 2 - 72ci„,C(3)) 
64 

+ Yg(-6c 2m + 15ci,„ + 24c 3m )x m + i(6c 2m - 3ci m )a;^ + ^c lm x^j e 2Xm 
+0( e ix ™) 

where Q is the Riemann Zeta function, and u m \ must match the outer solution 
which is given in the appendix by (|A.10[) and has only constant terms to this 
order. Hence we require c 2m = and c\ m — 0. The matched solution is now 

u mi( x m) = (1 - tanh 2 (a; m )) « m i (3-25) 
( _ 2e 6 ^ + 4 + lOe 2 - - 12e 4 ^ - 24x m e 2x A ^ , 

16 v J ( e 2x m + iy 

where we denote by u^ni the solution that is obtained by matching to the outer 
solution Y. As we will see later, exponential matching to the inner solution uo, 
i.e. as x m — > oo, where we find 



u; nl '(x m ) = -c 3m e m + -c 3m + I --c 3m + AK ml + -c 3m x m ) e 

^-c 3m - 8/Cmi - 3c 3m x m ^ e~ AXm + 0(e~ 6Xm ) , 

requires also c 3m = 0. Hence, denoting by the solution that has been 
exponentially matched to the inner solution uq near x — 0, we obtain 

u mi( x m) = (l - tanh 2 (x m )) k,„i . (3.26) 

Internal layer near x = The 0(S 1/3 ) proble m is 

£(u i,x ) = /oi(-^o) , (3.27a) 

uoi(0) = 0, <!(0) = and u^'(0)=0, (3.27b) 

with 

/oi(a;o) := cio^Q + c 20 a;o + c 30 . (3.28) 
Its general solution reads 

woi(xo) = ¥>oi(zo) -0o(a;o) , (3.29) 
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where the homogeneous solutions are as before and the inhomogeneous solution 
is given by equation (|3 . 2 1 [) with a = 0, j = 1 and 

7o = - j^cw + hi(2) c 20 - c 30 , (3.30) 

so that ifoi (0) = and ipoi grows algebraically as xq — > — oo. Furthermore, 
symmetry requires </?oi(0) = and (^"'(0) = 0, which implies c±q = and 
C30 = leading to 

m(«o) = i 6 (i + C2 p -2, o) 2 I 1 - 4a; o + 12dilog( e 2a: ° + l)e- 2a: ° - e- 4a ° + \2x\e^ 

+ ir 2 e- 2xo + 12x e- 4x ° - Uxoe~ 2x ° - ln(l + e - 2x °)e 2x ° + %e~ ixa ln(l + e~ 2x °) 
- 8 ln(l + e- 2xo )) + e- 6xo ln(l + e~ 2x °) + 2e- 6x °x ^ , (3.31) 

where dilog denotes the dilogarithm function. The remaining free parameters 
of uoi to be matched are C20 and g\. As will be demonstrated later, exponential 
matching to u m requires an expression for uoi as xo — * — 00 

uoi(xo) = -^-e~ 2x ° -\c2oXo-\gx (3.32) 
16 4 8 



+ —[2c 20 n z + 15c 20 + 26 5 i + (48 gi - I2c 20 )x + 24c 20 ^ J e 2x ° 

+ i( ~ 3651 ~ 89020 ~ 6C207r2 + (84C2 ° ~ 144 - 9l)x ° ~ 72c 20^o)e 4a;o + O{e 6xo ) 

and then re-expanding u in the variable a; m . This shows that also c 2 o = 0, 
<ji = and C3 m = 0. Any other choice leads to a system for the parameters 
having no solution. Hence, only n m remains as a free constant in the two regions. 
The exponentially matched solution is therefore simply 

u^o) =0. (3.33) 

3.1.3 0(S 2/3 ) 

Internal layer near n m The problem of order S 2 ^ 3 is 

) = fm2(x m ) i (3.34a) 

U m2 {0) = ~u' m0 (0) K m2 - -u'mOl<Ll - Mml(°)«ml = «m2 ~ u' ml (Q) K ml , 

(3.34b) 

where 

fm2(x m ) ■= di m x 2 m + d 2m x m + d 3m + 6 u m0 (ii„{) 2 ■ (3.35) 

Note that (0) = 0. Again we choose the inhomogeneous solution so that it 
grows only algebraically as x m — > —00 and vanishes at x m = to obtain (|3.21[) 
with a = to, j = 2 and 

7m2 = -j^di m + ln(2) d 2m - d 3m - n 2 ml , (3.36) 
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so that the general solution is represented as 

U m2 (x m ) = -K m 24> m {x m ) + tf m 2{x m ) . (3.37) 

As x m ^-cowe have to compare 

«m2(lm) = -^(dxm + 2 ^3m) ~ T^m^m ~ T^lm^m 
o 4 4 

+e 2 " m (^[(-7 - 72C(3))d lm - 8d 3m + 256(« m2 - K , 2 nl ) + (30 + An 2 )d 2m ] 

3 3 1 \ 

+ y^(5di m - 2d 2m + 8d 3m )x m + -(2d 2m - d lm )x 2 n + -d lm x 3 m J + 0(e iXm ) 

with the outer solution. Matching the constant and the linear terms in x m yields 

1111 23 
- = -Ax - -A 2 + -CUi + —C\ + D x , (3.38) 

- \d 2m = 2 1 '*C 1 . (3.39) 

There is no quadratic term in the outer solution (|A.10[) , hence d\ m = 0. There 
are further matching conditions but they do not simplify the problem struc- 
turally at this point and will be enforced later, so that d 2m , d 3m and K m2 
remain to be determined via exponential matching. As x m — > oo, the expansion 
to this order can be written as 

f 1 11 1 g-2a; m , 
U m2 = 2 d3m ~ 2 d2mXm + g d 3m e2Xm 32~ ( _ 5Qd,3m ~ 1M 2m (3.40) 

- 2d 2m ir 2 + I28(n 2 nl + n m2 ) + (48d 3m - \2d 2m )x m ~ 24d 2m x 2 n ") + 0{e~ Axm ) . 

Internal layer near x = As for the 0(<5 1 / 3 ) problem, at 0{5 2 ' z ) we have 

C{u 02l x ) = f 02 (xo) , (3.41a) 

u 02 (0) = 0, Uo' 2 (0) = and u£ 2 '(0)=0, (3.41b) 

with 

fo2(xo) ■= d w xl + d 2a x a + d 30 ■ (3.42) 

The general solution is 

uo 2 (x ) = <Po 2 (x ) + g 2 ip (x ) , (3.43) 

where the homogeneous component is as before and the inhomogeneous part 
is obtained by setting a = 0,j = 2 and 702 = in (|3.2ip . so that 9502(0) = 
and (/?o 2 grows algebraically as xq — > —00. Symmetry requires y ' 2 (0) = 0, 
tfiQ2 (0) = 0, which implies d±o = and d 3 o = 0. The remaining free parameters 
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to be matched are d 2 o and g 2 . In order to exponentially match to u rn to 0{5 2 ^) 
and obtain m^ 2 ; we again have to expand 1(02(^0) as xq — * — 00, giving 

"02(2:0) = -rre - -d 2a x a - -fi 
16 4 8 

+ — ((15 + 2tt 2 + 2 ln(2))d 2 o + 26.g 2 + (48/t - 12d 2 o)^o + 24d 20 Zo) e 



(3.44) 



+ — - (89 + 6^)d 20 - 36A + (84rf 20 - 144fi)x - 72d 20 x 2 )e ix ° + O(e bxo ) . 

and re-express in terms of x m , where we have used the abbreviation fi = 
d 2Q ln(2)+ 52 . 



3.1.4 0(S) 

Internal layer near K m The problem to be solved at order 0(6) is 
Um3(0) = -u' m2 (0)n m i - u mO (0)n m in m2 - u' m0 (0)K m3 

- g W mo(°) K ml - u ml(°)«m2 - ^<„i(0)^ n2 , 



with 



f m 3(x m ) ■= 2 ((u^) 3 + 6 U m0 
_ 2 3/2 



(3.45a) 
(3.45b) 

(3.46) 



-dilog(e 2a: ™ +!) + -(! + k lm )x 2 m + (ln(2) + k 2m )x m + k 3m 



Again we choose the inhomogeneous solution so that it grows only algebraically 

as Xm * 

and vanishes at x m = and so that we obtain (p m 3(x m ) by using 
formula (|3.21|1 with a = m 1 j = 3 and "f m 3 = 0. The solution is 



Um3(Xm) = ~ «m3 (0)(j) m (x m ) + lf m 3 (x n 



(3.47) 



where k\ m , k 2m , k 3m and n m3 remain to be determined via matching. In order 
to exclude exponential growth as x m -> -oo wc obtain the relation 



V2 



1 2 m 



(n,ni (-(12 + 9n 2 )d 2m + 12d 3m - 24K m2 ) 



81n(2) 

^V2(24fc 3m -121n(2) 2 + A :im 7r 2 ; 



(3.48) 



so that the expansion obtained 

1 



-oo IS 



u m 3{x m ) = J-^i 1 + + 4fc 3m) + -^=(ln(2) + k 2m )x Tl 
+ (k lm + l)?^x 2 m + 0(e 2 *-). 



(3.49) 
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Comparing this with the outer solution to 0(5), equation (|A.10p . yields the 
matching conditions 

^=(1 + k lm + 4fc 3m ) = \ -±A! + |ci) A 2 + (J^C 2 + ±D X ) A, (3.50) 

+ Ia 3 - —dA\ - — 2 1 /3C 1 + K x - —C l D 1 + —A\ + —C\ 
2 3 216 1 12 1 1 7 16 1 28 1 

for the constant terms, 
-L(hi(2)+fc 2m ) = (£> 1 -yC 1 2 )2 1 / 6 and (k lm + 1) ^ = 2" 2/3 Ci (3.51) 

for the linear and the quadratic terms, respectively. 
Expanding the solution as x m —* oo we find 

u m3 {x m ) = — {n m id 2m {^ 2 + 24) - 48K m id 3m + 2V27r 2 (l - k lm ) - 48V2k 3m j e 

+ — (n ml d 2m (Tl^ 2 + 72) + V2(fc lm (12 - 6tt 2 ) - 96fc 3m - 12 + 2ir 2 ) 
96 V 

+ -i=(ln(2) + k 2rn )x m + (k lm + l)^x 2 m + 0(e~ 2 ^) , (3.52) 

and wc will exponentially match it to the solution near x = 0, which we construct 
next. 

Internal layer near x = The general solution to the 0(5) problem 

£(uo3,xo) = /osfco) , (3.53a) 

u 03 (0) = 0, uo' 3 (0) = and u^'(0)=0, (3.53b) 

with 

/oaOro) := -I 1 ' 2 [dilog(e 2 * + 1) - dilog(2) + 2fi 2 x Q + (1 + k w )x 2 ] (3.54) 

and the abbreviation \i 2 = In (2) + ^20 reads 

"03(^0) = ^03(2:0) +g 3 ipo(x ) , (3.55) 

where we have required that uo 3 (0) = and w ' 3 (0) = 0. If we also enforce 
1*03(0) = then fcio = 0. Again we take an inhomogencous solution (^03(^0) 
which satisfies the above conditions, so that the general solution is obtained 
with 

f 1 1 ( z 2 + 1\ 2 lnf2z) 2 

Hi = v / 2(ln(2) 2 +2fc 20 ln(2))- 53 and u = - In — K -^-dz « 0.3094 , 

J z \ 2z J z 
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"03 



+ 1(36 M1 + V2(12 - n 2 )) + ^x + ^x 2 (3.56) 

1 



192 



156#i + V2[(19 - 24fc 20 )7r 2 - 15 - 288w - 180/i 2 ] 



' -24/ix + \/2(12/i 2 - 11)) x + (3 - 12 M2 ) x 2 , 



16 V ^ ' ""v-^ "V u ' 8 v ^ u yT 



e aiB0 +O(e* B0 ). 



For exponentially matching to u m this again has to be re-expressed in x m and 
combined with the corresponding expressions for woo, Uoi an d u 02 ■ This will 
be done in the next section. 

3.2 Exponential matching 

Now we have to match the rest of the solution u m (x m ) to the rest of the solu- 
tion uq(xq). This requires matching the exponential terms in addition to the 
algebraic terms, similarly to the procedure for the CCH equation, i.e. matching 
of the solution describing the internal layer near x = n m to the solution near 
x = requires expressing the variable xq in terms of x m (or vice versa). Recall 
again that xo — x m + R m /y2 and that R m < 0; the e 2x ° terms in the uq expan- 
sion will produce e 2Xm terms with a factor e^ 2 ~ Km (and analogously for e~ 2x ° 
terms) and so we will find their corresponding matching partner at a different 
order in 8 in the u m expansion, as we have shown for the CCH equation. The 
somewhat subtle difference here is that additionally we need to determine the 
relationship between e^ 2Km and 5 and we have in principle several choices, only 
one of which allows a consistent matching of both expansions. One can observe 
that the choice e^ 2Km = pS 1 ^ 3 , where p is some constant quickly leads to a 
contradiction. However, setting 

will lead to a 0(S 2 ^ 3 ) shift of terms, so that e.g. 

e 2xa will shift to a term 8 2/3 e 2Xm , (3.58) 

e' 2xo will end up as a term <T 2/3 e~ 2Xm (3.59) 

and so forth, so that e.g. a term e 2x ° in the leading order part of the uq expansion 
will have to match a e 2Xm term in the 0(S 2 ^ 3 ) part of the u m expansion, or a 
e~ 2x ° term in the 0(5) part of the uq expansion will have to match a e~ 2Xm 
term in the 0(S 1 ^ 3 ) part of the u m expansion. This will also produce terms 
that will have no partner term in the transformed expansion. Their coefficients 
must then be set to zero. If we now sum the expansions for uoi(xo), 1*02(^0) 
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and uo3(xo) and re-expand using (|3.57[) . we obtain 
!/..,_, \ 1 



u (x m ) 



-1 



16 



(d 20 ln(2)+. g2 )e- 2 ^p+ I ^( 



12/ii - V2tt 



2l' d 



48e^ 



^ ^36 M i + V2[12 + (161n(<5) - 241n(p))/x 2 + 6(ln(p) - | ln(<5)) 2 - tt 2 ]^) 



V2 \/2 
— (2 ln(<5) + 6^2 - 3 ln(p)) x m + —x 2 m 



<). 



which has to match u m i(x m ), u m 2(x m ) and u m 3(x m ) to each order, respectively. 
From this we obtain further conditions for the parameters in addition to those 
we have already found. Solving the complete system of equations then yields 
the solutions for the width of the hump 



A= V2, ( P_ 

6 \W(/3V 3 ) 3 



(3.61) 



with (3 = 2 11 / (27i5 2 ), where W is the Lambert W function (so W(x) is the solu- 
tion of x = W exp(W)). The expressions for the remaining matching constants 
C\,Di, etc. are omitted. The first correction in (|3.2p has the coefficient 



-3 2 1 / 6 . 



(3.62) 



Note that in the transformed expansions as well as in the expressions for the 
parameters also contain so-called logarithmic switch-back terms. 



(3.60) 



4 Numerical method for the fifth-order phase 
space 

For the numerical stationary solutions of the HCCH equation (|3.ip we apply 
the same scaling for u that we used for the CCH equation to obtain equilibrium 
points at ±1. 



(1-c 2 ) 



5VA 



lim c = =pl 

x — >-±oo 



(4.1) 



again assuming that derivatives vanish in the far field. Reduction to a first order 
system U' = F(U), with F : R 5 — > R 5 , gives a five-dimensional phase space, 
where the first four components of Fi(U) are equal to E/j+i and the fifth is 



F 5 (U) = 6A{U 2 y + I8AU1U2U3 + (3A(Uiy - l)C/ 4 + SVA(1 - {U x y)/2 . (4.2) 

The equilibrium points are U = ±(1, 0, 0, 0, 0) T and at these points the char- 
acteristic polynomials arc 



V ± (X) = X 5 + \ 3 {1-3A)±SVA. 



(4.3) 
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For small 5 the manifolds W U (U + ) and W S (U~) are both two-dimensional, re- 
sulting in a codimcnsion two event when searching for heteroclinic solutions 
connecting the two hyperbolic fixed points U + and U~ . The HCCH equation 
exhibits the same reversibility properties as its lower order version. This re- 
versibility is again given by the transformation (|2.8p from the CCH section, 
which also here fulfills RF(U) = —F(RU). The codimcnsion reduces by one 
and again we deal with a codimcnsion one problem and two parameters, hence 
we may expect solution branches in the {A, 5) parameter plane. Section 12.11 
showed that a condition for the existence of heteroclinic orbits is a value where 
the distance function (|2.9[) reaches zero and the same condition holds for the 
HCCH equation. The phase space is sketched in figure indicating the lin- 
earizations of the intersecting manifolds in the equilibrium points. For this 




Figure 6: HCCH: Heteroclinic orbits between the equilibrium points are sought 
in a 5-D phase space that is indicated here in 3D. The manifolds W U (U + ) 
and W S (U~) are two-dimensional which is suggested by the two planes in the 
picture. 

problem a shooting method will be very slow and may lead to bad accuracy 
since the additional parameter, say ip 6 [0, 27r), an angle defining points on a 
circle close to the equilibrium point on the linearization of the two-dimensional 
manifold, requires a very fine resolution to obtain heteroclinic solutions. 

4.1 Boundary value problem formulation 

There exist several possibilities to set up equations for finding heteroclinic con- 
nections in a boundary value problem framework. Generally one crucial stum- 
bling block is the choice of a suitable phase condition that picks a certain so- 
lution out of the infinitely many available ones due to phase shifts 0, Q . We 
choose to incorporate one phase condition proposed by Beyn 0], for which we 
use an approximation of the solution, V, typically given by a previous solution 
for slightly different parameter values. Equation (|4.1[) contains two parameters, 
A, 5, and in addition the truncated domain length L. As discussed by Doedel 
et al. [H one of the free parameters can be replaced by L to find a connection. 
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We replace 6, solve and continue after extrapolating to an approximate value of 
A for a nearby chosen and fixed 5. Rescaling the domain to [0, 1] yields, with 
the phase condition variable U p h introduced by Beyn Q the first order system 

U'i=LU i+u i= 1,2,3,4 (4.4a) 
U' 5 = L (6A(U 2 ) 3 + 18AUtU 2 U 3 + (SAiUr) 2 - l)C/ 4 + SVA ^ ~ jfi^ 



(4.4b) 

U' ph = Liy'fU (4.4c) 
L' = 0, A' = . (4.4d) 

Hence, we obtain one equation for the phase condition and two for the param- 
eters in addition to the five given by the original ODE, i.e., we have an overall 
system of eight equations which have to be supplemented by the same num- 
ber of boundary conditions. At the edges of the domain we utilize projected 
boundary conditions 0, Q , which make use of eigenvectors in the equilibrium 
points and can be incorporated by computing Vb, the matrix whose columns are 
composed by the eigenvectors which correspond to the eigenvalues at the upper 
equilibrium point U + with negative real part, and by forming the counterpart 
Vi containing those eigenvectors given by the unstable directions at the lower 
stationary point U~ . Hence, we consider the eight boundary conditions 

U ph (0)=0, U ph (l) = 0, V o T (U(0) -u+) = o, v 1 T (u(i)-u-) = o. 

(4.5) 

For initial estimates we can use solutions obtained from the asymptotic analysis 
of section I3.1[ i.e. the leading order solution tanh profiles 

V(x) = — tanh(a; — K) + tanh(a;) — tanh(:r + K) , 

for the het\ solution with guessed root-distance K. 

The boundary value solvers we use are based on mono-implicit Rungc-Kutta 
formulae [3, . As for the CCH problem efficiency can be improved by mak- 
ing use of the theory from section |2~T1 which holds analogously for the HCCH 
equation to obtain a boundary condition at the fixed point of a point-symmetric 
solution. We can use half of the previous domain length and phase conditions 
become redundant, because the phase is already fixed. We replace the projected 
boundary conditions by 

Ui{0) = l, U 2 (0f + U 3 (0) 2 -0, U 4 {0) 2 + U 5 (0) 2 = 

so that together with the self-reversibility condition on the right interval end 
Ui(l) = Us(l) = Uz,(l) = we have six conditions which match the five equa- 
tions together with the free parameter A. Final solutions are obtained by re- 
flecting the solution and its derivatives around zero and changing the signs of 
the first, third and fifth component. Examples of branches of different solutions 
are shown in figure El 
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Figure 7: (V~A, 5)— plane with curves for the first three heteroclinic connection 
branches for the HCCH equation. The dashed line in the parameter plane 
indicates the position where the positive roots of the characteristic polynomial 
in U + have nonzero imaginary parts. Below and to the right we see five phase 
space diagrams (tuples (Ui, U2), (Z/i, t/3), . . .) for selected solutions pointed out 
with arrows marking the corresponding parameters. The first pair (£7 1 , t/2 ) is 
plotted as bold solid curve. 



4.2 Solutions and comparison to analytical results 

With the boundary value formulation we are able to compute new HCCH sta- 
tionary solutions. In figure [8] we see a particular heti solution and the profile 
of the growing structure. 



profile 




Figure 8: heti solution for 5 = 0.01 and A = 0.443 and the corresponding profile 
obtained by integration. 

Up to 3D one can nicely visualize heteroclinic orbits in the corresponding 
phase space, while when the dimension is four or higher and the derivatives van- 
ish in the far field one can still plot the 2D phase spaces (Ux, U2), (£7i, U3), . . . 
and demand connections between the equilibrium tuples (±y/A, 0) as a neces- 
sary condition for heteroclinic orbits in the higher order space. Several such 
projections onto 2D are shown in figure [3 where we also see a very rapidly os- 
cillating heteroclinic curve in the bottom left plot which was found by a shooting 
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approach with a minimization procedure that used the two parameters and an 
angle as free parameters and the distance function (|2.9|) as objective function, 
depending on those parameters. It indicates that as shown for the CCH equa- 
tion we can in fact find many more hetk branches than those presented for 
k = 0, 1,2, all emerging from (A, 6) = (1,0), which corresponds to the Cahn- 
Hilliard equation. 

In figure[9]wc see the change in appearance of solutions on the het2 branch as 
S is increased. The shape varies from a solution with two pronounced humps to 
a monotone one, similar to the het$ solution, although associated with different, 
smaller, values of A. This is crucial if one wants to compute solutions for bigger 8 
with a boundary value solver. It easily happens that the solver switches between 
solution branches, however, this can be prevented by starting continuation in a 
parameter regime where the high-slope parts of the solutions are non-monotone, 
and continuing with small steps. A characteristic of the HCCH solutions is the 



8=0.002 8=0.005 8=0.015 8=0.02 




Figure 9: Structural change of the scaled het^ solution as 5 is increased. 

overshoot from the equilibrium value before the solutions go down. This is not 
observed for the CCH equation, where the shape is similar at these regions to 
hyperbolic tangent functions. 

In light of the expansion (|3.2p we try to estimate the 0(6 1 ^ 3 ) terms A\ 
for the different heteroclinic connections in a range of very small 6. As we 
see in figure [TU1 on the left, the numerically obtained values for A behave like 
A = 1 — 2 1 / 6 <5 1 / 3 in case of the het solutions, so that A\ = — 2 1 / 6 , which is 
consistent with the result in Savina et al. [22j. The numerical result for het\ is 
in line with the analytical value (|3.62|) and since for heti we see the agreement 
A\ ss —5 2 1 / 6 , we propose for higher order trajectories that for hetk we have the 
general approximation A\ ~ — (2k + 1) 2 1 / 6 , which is reminiscent of the CCH 
expression (]2 .45[> . Hence this formula is used in figure [TU] to plot the analytical 
values. 

We measure the distance between the first and second root for the het\ and 
the het2 solutions as seen in figure [10] on the right. We compare this to the 
analytical expression (|3.61|) for the one-hump solutions in the same figure and 
see that for small 6 the agreement is good. For both het\ and het^ solutions 
the distance is seen to increase logarithmically as 6 decreases. 
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Figure 10: Left figure: Logarithmic version of the (Vj4, S) plot for very small 
5. Drawn through curves giving the analytical values, dash-dotted lines those 
computed with the BVP solver. On the right we see the distances between 
the first two roots of the het\ , het2 solutions, numerically and for het\ via the 
analytical expression (|3 . 6 1 1) (solid line). 

5 Conclusion 

We have demonstrated that a sixth-order generalization of the convectivc Cahn- 
Hilliard equation admits multiple stationary solutions connecting constant val- 
ues. As for the fourth-order convective Cahn-Hilliard solution, these include a 
simple base solution, which is monotone for the CCH and "almost" monotone 
for the HCCH equation. More complex solutions, containing multiple humps, 
are also possible for each value of the forcing parameter 8, given particular 
values of the integration constant A(S). These non-monotone stationary solu- 
tions constitute an essential part of the solution structure for this higher-order 
Cahn-Hilliard type equation. We have demonstrated this via a numerical in- 
vestigation of the phase space in which we arc able to follow solution branches. 
For the simplest of the multi-humped solutions, the het\ branch, careful use of 
matched asymptotics that accounts for exponentially small terms allows us to 
find a solution which yields both the length scale for the solution (the "hump 
length") and the parameter value A(5) at which it occurs, in the limit of small 
d. Extension of the analysis to higher branches appears feasible. Our numerical 
evidence suggests that similarly simple asymptotic expressions hold for these 
branches, for both the CCH and HCCH equations. Physically, these solutions 
may represent situations where the edge energy regularization represented in 
(|1.4j) fails to produce a smooth transition between facets. 

Various issues, such as the stability of these solutions are presently being 
considered in the light of applications of the HCCH equation as a model for 
the morphology and dynamics of quantum dots. In particular, how do adjacent 
internal layers derived from these solutions interact, and what is their effect on 
the coarsening behavior in large spatial domains? Savina et al. [12] have begun 
an investigation of these questions by numerical simulation of (|1.4|) : it is likely 
that asymptotics can yield further insights. 

Physically, further interesting questions relate to the extension of the HCCH 
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model to richer models for the energetics of facetted surfaces, and analyzing the 
three-dimensional extension of the model. 
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A Outer Problem 

For the solution to the outer problem (|3.6|) . (|3.7p it is easy to observe that to 
leading order in 6 the solution of 

Q (Y ) = -V2 (Y 2 - 1) with lim F (0=1 (A.l) 

sss £ — ► — oo 

is 

y (O = i- (A.2) 

To 0(8 1 ^ 3 ) the general solution to the problem 

Y 1(U - V2F! = with ^^(0 = ^ (A.3) 



is 

A\ 2 1/6 



,5/6 



Yi(0 = ^ + Cie 217 ^ + e~«/ 25/6 [c 2 cos (^£/2 5 / 6 ) + C 3 sin (x/3 £/2 5 

(A.4) 

with Ci,C2 and C3 being constants of integration. The far field condition 
requires that Y\ remains bounded as £ — > —00. Hence, C2 = C3 = and 

n(£) = ^ + ^ 2l/6? (A.5) 

Using this and the far field conditions, the solution to the 0(<5 2 / 3 ) problem 

Y 2 ^-V2Y 2 = -^-±(3{Y 1 2 ) m -V2Y?) with ^Yatf) = 

(A.6) 

is 

tt(0 = t - f + ^ g2l/S ' + ^ ( x " 2l/6 ^) - l^ 2 "' (A7) 

and to the 0(<5) problem 

Y 3m -V2Y 3 = -^ + lW + V2y I K-l(l?+6Y I y J )« ( 

with ^y,®. + 4 (A.8) 
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it is 



Y 3 (0 



A 3 A X A 2 A\ 



4 16 

Ki - ^—C x + - (AiDi + AaCi) - i^C x A\ 



12 

'y/2„ 2 1 / 6 



216 



,/9 oi/B \ ol/3 

- v (^i^i + ^d) + ^2 1/6 ^d e + i- Afd e 



(A.9) 
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with another integration constant K±. Finally, wc obtain the asymptotic repre- 
sentation in terms of x m : 



Y(x m ) = 1 + 



d + -A, 



$1/3 



d 2 1 / 6 : 



- A? + - d Ai + Di - — d 2 + - A 2 

g 1 -r 3 1 1 -r i 14 1^ 2 2 



d 2 2V6 Xm + Dl 2 V* Xm + I d 2 1 /3 a $ i + f-i Ax + i C 



7 

7 1 \ 1 

_C? + - Di )A 1 + - A 3 
12 1 3 J 2 3 



2 \ 4 3 

i^C 1 A 2 --2 1 / 3 d 
216 1 12 



23 1 

- yd D 1 + -A\ 
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